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Silicene is a monolayer of silicon atoms forming a two-dimensional honeycomb lattice, which shares almost 
every remarkable property with graphene. The low energy structure of silicene is described by Dirac electrons 
with relatively large spin-orbit interactions due to its buckled structure. The key observation is that the band 
structure is controllable by applying the electric field to a silicene sheet. In particular, the gap closes at a certain 
critical electric field. Examining the band structure of a silicene nanoribbon, we demonstrate that a topological 
phase transition occurs from a topological insulator to a band insulator with the increase of the electric field. 
We also show that it is possible to generate helical zero modes anywhere in a silicene sheet by adjusting the 
electric field locally to this critical value. The region may act as a quantum wire or a quantum dot surrounded 
by topological and/or band insulators. We explicitly construct the wave functions for some simple geometries 
based on the low-energy effective Dirac theory. These results are applicable also to germanene, that is a two- 
dimensional honeycomb structure of germanium. 



I. INTRODUCTION 

Graphene, a monolayer honeycomb structure of carbon 
atoms, is one of the most important topics in condensed mat- 
ter physics 1 . One of the obstacles of graphene for electronic 
devices is that electrons can not be confined by applying ex- 
ternal electric field 2 . Thus, graphene nanostructures such as 
graphene nanoribbon 3 and nanodisk 4 have been considered, 
which are to be fabricated by cutting a graphene sheet. Re- 
cently a new material, a monolayer honeycomb structure of 
silicon called silicene, has been synthesized 5-7 and attracts 
much attention 8-10 . Silicene has Dirac cones akin to graphene. 
Almost every striking property of graphene could be trans- 
ferred to this innovative material. Furthermore, silicene has 
advantage of easily being incorporated into the silicon-based 
electronic technology. 

Silicene has a remarkable property graphene does not share: 
It is the buckled structure^ owing to a large ionic radius of 
silicon (FigQ]). Consequently, silicene has a relatively large 
spin-orbit (SO) gap of 1.55meV, as makes experimentally ac- 
cessible the Kane-Mele type quantum Spin Hall (QSH) ef- 
fect or topological insulator 9,10 . Topological insulator 11 -^ is 
a new state of quantum matter characterized by a full insu- 
lating gap in the bulk and gapless edges topologically pro- 
tected. These states are made possible due to the combination 
of the SO interaction and the time-reversal symmetry. The 
two-dimensional topological insulator is a QSH insulator with 
helical gapless edge modes 13 , which is a close cousin of the 
integer quantum Hall state. QSH insulator was proposed by 
Kane and Mele in graphene 14 . However, since the SO gap is 
rather weak in graphene, the QSH effect can occur in graphene 
only at unrealistically low temperature 1 ^^. 

The buckled structure implies an intriguing possibility that 
we can control the band structure by applying the electric field 
(FigQ}. In this paper, we analyze the band structure under the 
electric field E z applied perpendicular to a silicene sheet. Sil- 
icene is a Z2 topological insulator 9 &lE z =0. By increasing 
E z , we demonstrate the following. The gap decreases linearly 
to zero at a certain critical field E c and then increases linearly. 
Accordingly, silicene undergoes a topological phase transition 




FIG. 1: (Color online) Illustration of the buckled honycomb lattice 
of silicene. A honeycomb lattice is distorted due to a large ionic 
radius of a silicon atom and forms a buckled structure. The A and B 
sites form two sublattices separated by a perpendicular distance 2t. 
The structure generates a staggered sublattice potential in the electric 
field E z , which leads to various intriguing pheneomena. 



from a topological insulator to a band insulator. At the criti- 
cal point (E z = E c ), spins are perfectly spin-up (spin-down) 
polarized at the K (K') point. 

We also investigate the zero-energy states under an inho- 
mogeneous electric field E z (x,y) based on the low-energy 
effective Dirac theory. There emerge helical zero modes in 
the region where E z (x,y) = E c . It is intriguing that the 
region is not necessary one-dimensional: The region can be 
two-dimensional and have any shape, where Dirac electrons 
can be confined. It is surrounded by topological and/or band 
insulators. Our result may be the first example in which heli- 
cal zero modes appear in regions besides the edge of a topo- 
logical insulator. The region may act as a quantum wire or a 
quantum dot. We construct explicitly the wave functions de- 
scribing helical zero modes for regions having simple geome- 
tries. In conclusion, we are able to realize a dissipationless 
spin current anywhere in the bulk of a silicene sheet by tuning 
the electric field locally. 




FIG. 2: (Color online) The band gap A as a function of the electric 
field E z . The gap is open for E z / ±E C , where silicene is an 
insulator. It can be shown that it is a topological insulator for \E Z \ < 
E c and a band insulator \E Z I > E c . 



II. TOPOLOGICAL AND BAND INSULATORS 

Silicene consists of a honeycomb lattice of silicon atoms 
with two sublattices made of A sites and B sites. The states 
near the Fermi energy are tt orbitals residing near the K and 
K' points at opposite corners of the hexagonal Brillouin zone. 
We take a silicene sheet on the xy-plane, and apply the electric 
field E z (x,y) perpendicular to the plane. Due to the buckled 
structure the two sublattice planes are separated by a distance, 
which we denote by 2£ with i = 0.23 A , as illustrated in FigQ] 
It generates a staggered sublattice potential oc 2£E z (x, y) be- 
tween silicon atoms at A sites and B sites. 

The silicene system is described by the four-band second- 
nearest-neighbor tight binding model 10 , 
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The first term represents the usual nearest-neighbor hopping 
on the honeycomb lattice with the transfer energy t — 1.6eV, 
where the sum is taken over all pairs (i, j) of the nearest- 
neighboring sites, and the operator c\ a creates an electron 
with spin polarization a at site i. The second term rep- 
resents the effective SO coupling with Aso = 3.9meV, 
where a = (cr x , a y , a z ) is the Pauli matrix of spin, = 



di x dj 



with dj and d« the two nearest bonds 



(di x d^j j 

connecting the next-nearest neighbors, and the sum is taken 
over all pairs of the second-nearest-neighboring sites. 

The third term represents the Rashba SO coupling with Ar = 
0.7meV, where fiij = ±1 for the A (B) site, and = 

dij j dij . The forth term is the staggered sublattice potential 

term, where Q = ±1 for the A (B) site. Note that the first and 
the second terms constitute the Kane-Mele model proposed to 
demonstrate the QSH effect in graphena^. 

The same Hamiltonian as (I2.lt can be used to describe 
germanene, that is a honeycomb structure of germanium^, 
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FIG. 3: (Color online) Band structure of silicene at the critical elec- 
tric field E c . (a) A bird's-eye view. Dirac cones are found at 6 cor- 
ners of the hexagonal Brillouin zone, (b) The cross section contain- 
ing a pair of K and K' points. The solid red (dashed blue) band is for 
up-spin (down-spin) electrons, which are gapless (gapped) at the K 
point but gapped (gapless) at the K' point. 



where various parameters are t = 1.3eV, Aso = 43meV, 
Ar = 10.7meV and t = 0.33 A. Hence the following anal- 
ysis is applicable to germanene as well. 

We study the band structure of silicene by applying a uni- 
form electric field E z . By diagonalizing the Hamiltonian 
(I2.lt , the band gap A (E z ) is determined to be 



A(E z ) = 2\£E z - V s z X : 



sol 



(2.2) 



where s z = ±1 is the electron spin and rj = ±1 is for the 
K or K' point (to which we refer also as the K± point). See 
also the dispersion relation (13.4b which we derive based on the 
low-energy effective theory. We emphasize that it is indepen- 
dent of the Rashba SO coupling Ar. The gap (12.21) vanishes at 
E z = V s zE c with 



E c = \so/t = 17meV/A. 



(2.3) 



We plot the band gap A (E z ) in the Fig J2] 

The gap closes at E z = ±E C , where it is a semimetal due 
to gapless modes. We show the band structure at E z = E c in 
FigOl It follows from (12.21) that up-spin (s z = +1) electrons 
are gapless at the K point (rj = +1), while down-spin (s z = 
— 1) electrons are gapless at the K' point (rj = —1). Namely, 
spins are perfectly up (down) polarized at the K (K') point 
under the uniform electric field E z = E c . 

It follows from the gap formula (12.2b that silicene is an in- 
sulator for E z 7^ ±E C . In order to tell the difference between 
the two insulators realized for \E Z \ < E c and \E Z \ > E c ,we 
study the band structure of a silicene nanoribbon with zigzag 
edges. The gap structure is depicted at two typical points, 
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FIG. 4: (Color online) One-dimensional energy bands for a silicene 
nanoribbon. (a) The bands crossing the gap are edge states, demon- 
strating that it is a topological insulator. There are two edge states 
since a nanoribbon has two edges (red and blue lines for the left and 
right edges), (b) All states are gapped, demonstrating that it is a band 
insulator. 



E z = and E z = 2E c ,'m Fig 01 We see that there are gap- 
less modes coming from the two edges at \E Z \ < E c , as is the 
demonstration of a topological insulator-. On the other hand, 
there are no gapless edge modes for \E Z \ > E c , showing that 
it is a band insulator. We conclude that a topological phase 
transition occurs between a topological insulator (\E Z | < E c ) 
and a band insulator (\E Z \ > E c ) as E z changes. 

The reason why gapless modes appear in the edge of a topo- 
logical insulator is understood as follows. The topological in- 
sulator has a nontrivial topological number, the Z2 index 14 , 
which is defined only for a gapped state. When a topological 
insulator has an edge beyond which the region has the triv- 
ial Z2 index, the band must close and yield gapless modes in 
the interface. Otherwise the Z2 index cannot change its value 
across the interface. 



III. LOW-ENERGY DIRAC THEORY 

We proceed to analyze the physics of electrons near the 
Fermi energy more in details. The low-energy Dirac theory 
has been proved to be essential in the study of graphene 17 and 
its various derivatives^*^. It must also be indispensable to ex- 
plore deeper physics of helical zero modes and promote fur- 
ther researches in silicene. 

We may derive the low-energy effective Hamiltonian from 
the tight binding model (12.1b around the K v point as 10 



Hrj = Hv F (k x r x - rjkyTy) + rjTzhu + £E z r z , (3.1) 



hn = -Xso^z - «A R (k y (j x - k x (j y ) , 



with 

(3.2) 

where r a is the Pauli matrix of the sublattice, = ^-at = 
5.5 x 10 5 m/s is the Fermi velocity, and a = 3. 86 A is the lattice 
constant. It is instructive to write down the Hamiltonian H+ 
explicitly as 
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in the basis {^At> ^Bt> i>Ai, where k± = k x ±ik y . 

The two Hamiltonians H+ and H- are related through the 
time-reversal operation. 

The energy spectrum is readily derived from (13. IK 



£ v = ±\Jh 2 v$k 2 + [£E Z - rjs z ^ Af + a 2 A|/c 2 J , (3.4) 
which yields the result (12.21) to the gap energy at k = 0. 

IV. INHOMOGENEOUS ELECTRIC FIELD 

The low-energy Dirac theory allows us to investigate ana- 
lytically the properties of the helical zero mode under inhomo- 
geneous electric field. In so doing we set Ar = to simplify 
calculations. This approximation is justified by the following 
reasons. First all all, we have numerically checked that the 
band structure is rather insensitive to A R based on the tight- 
binding Hamiltonian d2.lt . Second, A R appears only in the 
combination X^k± in the Hamiltonian d3.lt , which vanishes 
exactly at the K± points. Third, the critical electric field E c is 
independent of Ar as in (12.3b . 

A. Inhomogeneous electric field along x-axis 

We apply the electric field E z (x) perpendicularly to a sil- 
icene sheet homogeneously in the y direction and inhomoge- 
neously in the x direction. We may set k y =constant due to 
the translational invariance along the y axis. The momentum 
k y is a good quantum number. Setting 



^ (x,y) = e ikyV $ (x) . 



(4.1) 



we seek the zero-energy solution, where ^ (x, y) is a four- 
component amplitude. The particle-hole symmetry guaran- 
tees the existence of zero-energy solutions satisfying the rela- 
tion (j>B (x) = i£(j)A (x) with £ = ±1. Here, <p A is a two- 
component amplitude with the up spin and the down spin. 
Then the eigenvalue problem yields 



H v (t) A (x) = E^(j) A (x), 
together with a linear dispersion relation 

Er)£ = T)£hv F k y . 

The equation of motion for (j) A (x) reads 

(£Hv F d x + Tf\so(T z - IE Z (x)) <j>A{x) 
We can explicitly solve this as 

(j) ASz (x) = Cexpf(x), 

with 



0. 



(4.2) 



(4.3) 



(4.4) 



(4.5) 



f(x) = Cexp 



(3.3) 



(-ris z Xso+£E z (x'))dx f 



(4.6) 
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where C is the normalization constant. The sign £ is de- 
termined so as to make the wave function finite in the limit 
\x\ — ^ oo. The current is calculated as 



J Sz (x) Re 



2mi Sz y 6 



^ Asz (x)\\ (4.7) 
m 



This is a reminiscence of the Jackiw-Rebbi mode 20 proposed 
for the chiral mode. 

The difference between the chiral and helical modes is the 
presence of the spin factor s z in the wave function. As we 
shall see explicitly in some examples in what follows, we find 
the condition either rjs z = 1 or rjs z = —1 for convergence 
of the wave function. The condition r]s z = 1 implies that the 
spin is up (s z = 1) at the K point (rj = 1) and that the spin is 
down (s z — —1) at the K' point (rj = — 1). Consequently, the 
up-spin electrons flow into the positive x-direction while the 
down-spin electrons flow into the negative x-direction, im- 
plying that the pure spin current flows into the positive in- 
direction. On the other hand, the condition r\s z = — 1 implies 
that the pure spin current flows into the negative ^-direction. 

Interface between topological and band insulators: We ap- 
ply an electric field such that 



E z (x) = olxIL 
Substituting it to (14.6b . we obtain 



\x) 



2hv F 



r)s z \ 



so 



r]s z \$o 



a 



(4.8) 



(4.9) 



where £ is chosen to make £a < for convergence. The wave 
function is localized along the two lines x = r]s z Xso/ct = 
±\so/ot taken in the bulk, where the gapless mode emerges 
since E z (x) = ±Aso- The currents are helical along these 
two lines, where r]s z = 1 along one line and r]s z = — 1 along 
the other line: The spin currents flow in the opposite directions 
along these two lines, sandwiched by a band insulator (\x\ > 
Aso/M) and a topological insulator (\x\ < \so/\a\). We 
illustrate the probability density \(j)As z (%) | 2 in FigOa) in the 
case a > 0. Each line may be used as a quantum wire. 

Interface between metal and insulator: We apply an electric 
field such that 



ax , 



E z ( x )=E c + —®{x), 



(4.10) 



where 6 (x) is the step function: 9 (x) = for x < and 
6 (x) = 1 for x > 0. Substituting it to (14.6b . we obtain 



/(*) 



hw \ 2 v 



) + (1 - rjs z ) \ so x 



Y 



(4.11) 



where £ is chosen to make £a < for convergence. Namely, 
when we choose a > 0, we find 



E z (x) = E c (metal) for x < 
E z (x) > E c (BI) for x > 



(4.12) 
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FIG. 5: (Color online) The probability density of the helical zero 
mode under the electric field E z . (a) E z is given by (14.81 ) with a > 0. 
The spin current flows between a topological insulator and a band 
insulator, (b) E z is given by ( 14.101 ) with a > 0. The spin current 
flows in the metallic region. The arrows indicate the helical mode 
with spin up at the K point and with spin down at the K' point. 



and when we choose a < 0, we find 

E z (x) = E c (metal) for x < 
E z (x) < E c (Tl) for x > 



(4.13) 



where Tl and BI stand for topological and band insulators. In 
the both cases the region x < is a metal. Furthermore, it is 
necessary that r]s z — 1 for convergence of the wave function, 
as implies that the current is helical in the metallic region. We 
illustrate the probability density \<j>A\ (%) 1 2 in FigHfb) in the 
case a > 0, which is a constant for x < 0. 

It is intriguing that there emerge helical zero modes in 
metal. This is not surprising because the edge of a topological 
insulator is a sufficient condition but not a necessary condi- 
tion for the emergence of helical zero modes. The helical zero 
mode requires the massless Dirac fermion, the time-reversal 
symmetry and the spin-orbit interaction. Our example shows 
explicitly that it can appear in regions besides the edge of a 
topological insulator provided these conditions are satisfied. 



B. Inhomogeneous electric field along r-axis 

We apply a cylindrical symmetric inhomogeneous electric 
field E z (r) to a silicene sheet. The equation reads 

-Aso^, + tE z (r) hvpe" (id r - ±d G ) \{*I>a\ =[) 
hv ¥ e- iG (id r + ±d e ) Xsocr z - IE Z (r) ) \ ) 

(4.14) 

We solve this for zero-energy states by setting 



(4.15) 



Mr,0) J \e-^ e ^ B (r) J' 



The equation of motion Hxip = is transformed into 



fry* ( d r + — ) + - £E Z (r) ) <f> A = 0, (4.16) 
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FIG. 6: (Color online) The probability density of the helical zero 
mode under the electric field E z with a > 0. (a) E z is given by 
( |4.19l >. The spin current flows between a topological insulator and a 
band insulator by encircling the topological insulator, (b) E z is given 
by ( 14.211 ) with a > 0. The spin current flows in the disklike metallic 
region confined within a topological or band insulator. The arrows 
indicate the helical mode with spin up at the K point and with spin 
down at the K' point. 



which we solve as 



C 



0a Sz (r) -=exp/(r), 
\/r 



(4.17) 



with 



) = ^f(- 



-T]s z Xso+eE z (r'))dr', (4.18) 



where C is the normalization constant and £ = ±1. The sign 
£ is determined so as to make the wave function finite in the 
limit r — » oo. 

Interface between topological and band insulators: We ap- 
ply an electric field such that 



E z (r) = ar/L 
Substituting it to (14.18b . we have 



(4.19) 




where £ is chosen to make £a < for convergence. The wave 
function is localized along the circle r = f]s z \so/ct > 0, 
where E z (r) = ±Aso- When we choose a > it is necessary 
that r]s z = 1, and when we choose a < it is necessary that 
r]s z = —1. In any of the two cases, there emerges helical 
zero modes and the spin current flows along the circle between 
a topological insulator (r < Aso/M) and a band insulator 
(T > Aso/M). The direction of the spin current is opposite 
for a > and a < 0. We illustrate the probability density 



|^As z (r) | 2 in FiglSta) in the case a > 0. This region may be 
used as a quantum wire. 

Interface between metal and insulator: We apply an electric 
field such that 



Ez(r) 



a(r - R) 



9 (r - R) , 



(4.21) 



where © (r) is the step function. Substituting it to (14.181 ). we 
have 



f(r) 



-U 



2£ 



Q(r-R) + (l- V 8 z )Xso{r-R) 



(4.22) 

where £ is chosen to make £a < for convergence. It is 
notable that E z (r) = E c for r < R and hence the system 
is metallic there. The wave function describes an interface 
between a metal for r < R and an insulator for r > R. The 
insulator is a topological insulator when we choose a < and 
a band insulator when we choose a > 0. Since it is necessary 
that i]s z = 1 for convergence of the wave function, the current 
is helical in the metallic region (r < R). We show the proba- 
bility density \(j>A\ (0 1 2 in FigIMb), where cj>A-\ (r) =constant 
for r < R. This region may act as a quantum dot. 



V. CONCLUSIONS 

Taking advantage of the buckled structure of silicene we 
have demonstrated that we can control its band structure by 
applying the electric field E z . Silicene undergoes a topo- 
logical phase transition between a topological insulator and 
a band insulator as E z crosses the critical point ±E C . It is a 
semimetal at E z = ±E C . 

A novel phenomenon appears when we apply an inhomo- 
geneous electric field. We have explicitly constructed wave 
functions of helical zero modes for simple geometrical regions 
based on the low-energy effective Dirac theory. The results 
imply in general that helical zero modes can be confined in 
any regions by tuning the external electric field locally to the 
critical field (12.3b , E z (x,y) = E c . Our system may be the 
first example in which helical zero modes appear in regions 
besides the edge of a topological insulator. It is to be em- 
phasized that we can apply an inhomogeneous electric field 
so that a single silicene sheet contains several regions which 
are topological insulators, band insulators and metals. Such a 
structure may open a way for future spintronics. Our results 
are also applicable to germanene, that is a two-dimensional 
honeycomb structure made of germanium. 
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